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Abstract 

We  present  an  improvement,  by  nearly  an  order  of  magnitude,  of  a  recent  algorithm 
by  Papadimitriou  and  Silverberg  for  optimal  motion  planning  for  a  ladder  [PS].  The 
improved  algorithm  runs  in  time  0(n3a(n)  log2  n),  where  n  is  the  number  of  obstacle 
edges  and  where  q(??)  is  the  extremely  slowly  growing  inverse  Ackermann  function. 


In  [PS],  Papadimitriou  and  Silverberg  consider  the  following  restricted  version  of 
optimal  motion  planning.  Suppose  we  are  given  a  line  segment  B  =  PQ  (a  "ladder") 
free  to  move  in  a  two-dimensional  closed  polygonal  region  V  bounded  by  a  total  of 
77  edges.  Each  placement  of  B  can  be  specified  by  (x,6),  where  x  is  the  position  of 
P  and  9  is  the  orientation  of  PQ.  A  placement  (x,9)  is  free  if  at  this  placement  B  is 
fully  contained  in  V.  Let  (s,0a),  {t,0t)  be  two  given  initial  and  final  free  placements 
of  B.  The  goal  is  to  plan  a  collision-free  motion  of  B  from  {s,6s)  to  (t,0t)  having  the 
following  properties: 

(i)  The  endpoint  P  moves  along  straight  segments  connecting  pairs  of  corners  of  V 
and  not  intersecting  Vc  (such  pairs  are  said  to  be  visible  from  one  another:  s 
and  t  are  also  considered  as  corners  in  this  context). 

(ii)  The  total  distance  traveled  by  P  is  minimum  among  all  collision-free  motions 
from  (s,0s)  to  (i,0t)  satisfying  (i). 

The  algorithm  presented  in  [PS]  runs  in  time  0(n4logn).  It  proceeds  roughly 
as  follows.  It  first  computes  the  visibility  graph  VG  of  the  given  obstacle  corners 
(including  5  and  i)  in  time  0(n2)  (the  graph  connects  by  edges  all  pairs  of  corners 
that  are  visible  from  one  another).  Next,  for  each  edge  e  =  uv  of  VG.  it  spends 
0(n2  log??)  time  to  compute  all  possible  free  motions  in  which  P  moves  along  e  from 
u  to  v,  and  associates  an  appropriate  cost  with  each  such  motion.  Combining  the 
transitions  along  each  edge  of  VG  yields  an  overall  '"connectivity  graph"  and  the 
desired  optimal  motion  can  then  be  obtained  by  applying  Dijkstra's  algorithm  to  the 
connectivity  graph. 

In  this  note  we  show  that  the  step  of  processing  a  single  edge  of  VG  can  be 
performed  in  only  0(na(??)  log2  7?),  where  a{n)  is  the  extremely  slowly-growing  func- 
tional inverse  of  Ackermann's  function.  This  leads  to  an  overall  0(n3a(n)  log2  n) 
algorithm,  improving  the  result  of  [PS]  by  nearly  an  order  of  magnitude.  To  obtain 
this  improvement  we  use  tools  that  have  recently  been  developed  in  [GSS],  [EGPPSS] 
and  [Si]  for  the  analysis  and  calculation  of  a  single  face  (or  certain  collections  of 
faces)  in  arrangements  of  curves  in  the  plane.  One  of  the  purposes  of  this  note  is 
to  demonstrate  the  significance  and  broad  applicability  of  these  tools,  although  the 
application  at  hand  is  not  completely  routine.  We  suspect  that  a  scan  through  the 
literature  would  yield  many  other  problems  whose  current  solutions  can  be  improved 
by  using  this  new  machinery. 

To  explain  how  this  improvement  is  obtained,  we  first  review  in  more  detail  the 
machinery  used  in  [PS]  to  analyze  and  solve  the  subproblem  of  moving  P  along  an 
edge  e  =  uv  of  the  visibility  graph.    This  restricted  motion  has  onl}-  two  degrees  of 


freedom,  and  any  (free  or  non-free)  placement  of  B  with  P  on  e  can  be  specified  by 
the  two  parameters  {x,6),  where  x  is  the  length  of  uP  and  9  is  the  orientation  of  B 
(for  simplicity,  and  without  loss  of  generality,  let  us  assume  that  e  lies  on  the  .r-axis). 
The  parametric  space  of  all  these  placements  is  a  cylinder  Ce  =  {(x,0)  :  0  <  x  < 
L.  0  <  0  <  2tt},  where  L  =  \uv\. 

Each  obstacle  edge  or  corner  o  induces  a  contact  curve  -y0  in  Ce  that  is  the  locus 
of  all  placements  (x,6)  of  B  in  Ce  at  which  its  endpoint  Q  touches  o,  if  o  is  an  edge, 
or  at  which  B  passes  through  o.  if  o  is  a  corner.  Let  A  denote  the  arrangement  of 
the  contact  curves  in  Ce.  This  is  the  cylindrical  map  obtained  by  drawing  all  these 
curves;  its  vertices  are  endpoints  or  intersection  points  of  contact  curves,  its  edges 
are  maximal  connected  portions  of  the  contact  curves  not  containing  a  vertex,  and 
its  faces  are  maximal  connected  cells  not  meeting  any  contact  curve.  Any  motion  of 
B  within  Ce  which  starts  at  some  free  placement  z  and  remains  free  cannot  cross  any 
of  the  contact  curves  ~)0  and  thus  must  remain  within  the  face  /  of  A  containing  z. 
Conversely,  any  placement  in  /  is  clearly  reachable  from  z  by  a  collision-free  motion 
of  B  within  Ce.  See  [GSS]  for  more  details. 

To  apply  this  crucial  observation,  consider  all  free  placements  of  B  with  P  touching 
u  (within  Ce  this  is  the  circle  x  —  0).  As  noted  in  [PS],  the  orientations  of  these 
placements  fall  into  0(?r)  circular  arcs,  each  of  which  is  delimited  on  both  sides  by  a 
placement  that  is  an  intersection  of  some  contact  curve  with  x  —  0.  The  algorithmic 
step  that  concerns  us  is  to  determine,  for  each  one  of  these  free  arcs,  6,  which  free 
arcs  at  x  =  L  are  reachable  within  Ce  from  8  (more  precisely,  from  {0}  x  8),  and 
what  is  the  cost  of  reaching  each  of  them.  (Recall  that  this  cost  is  the  total  distance 
traveled  by  P;  if  P  does  not  backtrack  during  the  motion,  the  cost  is  simply  I,  but 
when  backtracking  is  allowed  the  cost  can  be  higher.) 

The  preceding  observations  therefore  suggest  the  following  approach. 

(a)  Compute  all  the  faces  of  A  that  meet  the  circle  x  =  0  (this  collection  is  known 

as  the  zone  of  x  =  0  [Ed].  [EGPPSS]). 

(b)  Decompose  each  such  face  /  into  .r-monotone  regions  by  drawing,  from  each  point 

on  df  that  is  locally  x-extremal,  vertical  segments  that  are  extended  within  / 
until  they  meet  again  df. 

(c)  Build  a  connectivity  graph  whose  nodes  are  all  these  regions  and  whose  edges 

connect  pairs  of  adjacent  regions  (along  a  common  vertical  segment).  Associate 
with  each  region  a  cost  equal  to  its  x-width. 

(d)  Combine  all  the  connectivity  graphs,  over  all  edges  e  of  VG,  into  a  single  overall 

connectivity  graph  CG.  Then  search  for  the  minimum  weight  path  within  CG 


between  the  initial  and  final  placements  of  B. 

The  first  step  can  be  accomplished  by  adapting  the  algorithm  of  [GSS]  (see  also 
[EGPPSS])  for  the  calculation  of  a  single  face  in  an  arrangement  of  curves.  In  step 
(a)  we  need  to  calculate  more  than  a  single  face,  but  all  the  desired  faces  belong  to 
a  single  zone,  so  that,  using  the  techniques  of  [EGPPSS],  we  can  regard  all  of  them 
as  a  single  face.  (Intuitively,  we  regard  Ce  as  an  infinite  cylinder,  but  clip  the  curves 
%  so  that  they  do  not  extend  to  x  <  0;  this  makes  all  the  desired  faces  portions  of 
a  common  single  unbounded  face.  Note  also  that  the  analysis  of  [GSS],  [EGPPSS] 
is  normally  carried  out  for  planar  arrangements;  however,  as  discussed  in  [GSS],  it  is 
easy  to  cut  our  cylindrical  arrangement  into  two  planar  arrangements,  carry  out  the 
analysis  in  each  planar  patch  separately,  and  then  glue  together  the  resulting  faces 
to  obtain  the  desired  final  faces.) 

The  analysis  of  [GSS].  [EGPPSS]  shows  that  the  total  combinatorial  complexity 
of  all  the  faces  in  question  is  0(Xs+2(nj),  where  s  is  the  maximum  number  of  inter- 
sections between  any  pair  of  contact  curves  (we  can  actually  use  a  slightly  smaller  s 
-  see  below),  and  where  A5+2(n)  is  the  maximum  length  of  (n,s)  Davenport-Schinzel 
sequences  [HS],  [ASS],  which  is  almost  linear  in  n  for  any  fixed  s.  Moreover,  it  is 
shown  in  [GSS]  that  all  these  faces  can  be  computed  in  total  time  0(As+2(n) log   n). 

Before  continuing  with  the  analysis  of  the  remaining  steps,  let  us  analyze  the 
value  of  5.  In  a  naive  approach,  the  maximum  number  of  intersections  between  any 
pair  of  contact  curves  is  easily  seen  to  be  2.  However,  we  can  reduce  this  number 
to  1  by  applying  a  recent  observation  of  [Si].  Specifically,  we  partition  the  collection 
of  contact  curves  into  three  subcollections.  The  first  subcollection  Ti  consists  of  all 
curves  ~y0  where  o  is  a  corner.  The  second  subcollection  r2  consists  of  the  portions 
of  all  contact  curves  -y0  where  o  is  an  edge,  which  lie  in  the  range  —  |  <  6  <  |,  and 
the  third  subcollection  T3  consists  of  the  portions  of  the  latter  curves  which  lie  is  the 
range  |  <  9  <  y-.  It  is  now  straightforward  to  verify  the  following 

Claim  1    Any  pair  of  curves  in  the  same  subcollection  Tj  intersect  at  most  once. 

This  claim  implies  that  if  we  form  the  arrangement  Aj  of  the  curves  in  Tj  then  the 
complexity  of  the  zone  of  x  =  0  in  Aj  is  only  0(A3(n))  =  0(na(n))  [HS].  and  it  can 
be  calculated  in  time  0(na(n)  log2  n).  Then  the  combination  lemma  of  [GSS]  implies 
that  if  we  superimpose  the  three  arrangements  Aj  on  each  other,  the  complexity  of 
the  resulting  zone  of  x  —  0  remains  0{nc\{r\\).  and  it  can  be  computed  from  the  three 
corresponding  zones  in  the  -4/s  in  additional  time  0(na(n)  log  v). 

Continuing  the  analysis  of  the  algorithm,  we  note  that  step  (b)  can  be  done  in 
time  0(??q(?7)  log  v),  using  a  simple  line  sweeping  technique  (as  a  matter  of  fact,  the 


algorithm  in  [GSS]  produces  this  vertical  decomposition  for  free).  Step  (c)  is  also 
trivial  to  accomplish,  in  time  0(na{n)).  Thus  the  complexity  of  the  overall  graph 
CG  is  0(n3a(n)),  and  it  can  be  searched  within  the  same  time  bound.  Thus  the 
overall  cost  of  step  (a)  dominates  all  the  other  steps,  and  we  can  conclude 

Theorem  1  One  can  plan  optimal  restricted  motion  of  a  ladder  as  above  in  time 
0(n3a(n)  log2  n). 

We  conclude  with  a  few  open  problems.  First,  [PS]  also  considers  optimal  motion 
planning  where  the  cost  of  a  path  is  a  linear  combination  of  the  distance  traveled 
by  P  and  the  total  amount  of  rotation  executed  by  D.  For  this  cost  measure.  [PS] 
presents  an  0(n5)  algorithm,  and  it  would  be  nice  to  improve  this  complexity  using  our 
techniques.  A  second  problem  is  to  extend  our  techniques  to  the  case  (not  studied  in 
[PS])  where  the  motion  of  P  is  unrestricted  (and  the  cost  is  just  the  distance  traveled 

by  P). 
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